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Multifractal analysis of the metal-insulator transition in the 3D Anderson model II: 

Symmetry relation under ensemble averaging 

Alberto Rodriguez,^- ^ Louella J. Vasquez,^ and Rudolf A. R6mer^'[j 

^Department of Physics and Centre for Scientific Computing, University of Warwick, Coventry, CV4 7AL, United Kingdom 
^ Departamento de Fisica Fundamental, Universidad de Salamanca, 37008 Salamanca, Spain 

We study the multifractal analysis (MFA) of electronic wavefunctions at the localisation- 
delocalisation transition in the 3D Anderson model for very large system sizes up to 240 . The 
singularity spectrum /(a) is numerically obtained using the ensemble average of the scaling law for 
the generalized inverse participation ratios Pq, employing box-size and system-size scaling. The va- 
lidity of a recently reported symmetry law [Phys. Rev. Lett. 97, 046803 (2006)] for the multifractal 
spectrum is carefully analysed at the metal-insulator transition (MIT). The results are compared 
to those obtained using different approaches, in particular the typical average of the scaling law. 
System-size scaling with ensemble average appears as the most adequate method to carry out the 
numerical MFA. Some conjectures about the true shape of /(a) in the thermodynamic limit are also 
made. 

PACS numbers: 71.30.-|-h,72.15.Rn,05.45.Df 
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I. INTRODUCTION 

The multifractal analysis (MFA) of electronic wave- 
functions IV'iP at the localisation-dclocalisation transi- 
tion has been a subject of intense numerical study for the 
last twenty yearsi^i^iii^i^iii^i^iiaiiiiHaadij^ r^.^^ ^p^ ^f 

the critical jV'iP in a system with volume L'^ is based on 
the scaling of the generalized inverse participation ratios 
(gIPR), Pq{l) = SfcMfc(0' defined from the integrated 
measure fJ.k{l) = J2i IV^^P ™ ^^1 ^i boxes with linear 
size I covering the system. At criticality the scaling law 
Pq(A) ex A'^(') is expected to hold in a certain range of 
values for A = 1/ L. The well-known singularity spectrum 
/(a) is defined from the T{q) exponents via a Legendre 
transformation, /(a^) = qoiq — T(q) and a^ = t' [q). The 
physical meaning of the /(a) is as follows. It is the frac- 
tal dimension of the set of points where the wavefunction 
intensity obeys \ipi\'^ ~ -^~"j that is in a discrete system 
the number Na of such points scales as L^^"\ 

Recently, the report of a remarkable analytical result 
concerning the existence of an exact symmetry relation 
in the f[a)^ has required a profound revision of all the 
techniques involved in the numerical MFA. The reported 
symmetry lawii^ requires Aq = Ai_g in terms of the 
anomalous scaling exponents, which can be obtained by 
Aq = r(g) — d{q — 1). The symmetry can also be written 
as 

aq + ai^q = 2d, (1) 

or in terms of the singularity spectrum itself as 



5f{a) = |/(2d - a) - [/(a) + d - «]| = 0. 



(2) 



At the metal-insulator transition (MIT) the f{a) is a 
convex function of a ^ 0, with a maximum at ao ^ d 
where f{ao) = d. The values of /(a) are not restricted 
to be positive.— li^ The symmetry ([2]) implies that the 
upper bound is a ^ 2d and that the values of f{a) for 
a < d can be mapped to the values for a > d, and vice 




FIG. 1; Pictorial representation of the general features of the 
multifractal spectrum at criticality. The dotted black areas 
highlight forbidden regions for f{a). The grey shaded area 
can be connected to the white area through the symmetry 
relation ([1},(I2]), and vice versa. These two regions are de- 
termined by different wavefunction amplitudes, |'0ip > L~'' 
(white) and IV'iP < ^ "* (light grey). The properties of the 
spectrum at the points corresponding to g = 1/2, and the 
symmetry-related q — 0, and g = 1 are explicitly included: 
f{ao) = d, f{ai) = Qi, /'(ai) = 1, fiai/2 = d) = d ~ Ai/a- 



versa. A pictorial sketch of the expected properties of 
/(a) is shown in Fig.[TJ 

So far the symmetry has been numerically found in dif- 
ferent critical models below 3 D^^i^^ . In a previous work 
the role of the symmetry law in the MIT for the 3D An- 
derson model was thoroughly studied by the authors us- 
ing the typical average of the scaling law for the glPRi^ 
This has been usually regarded as the preferred way to 
perform MFA. In this work we study the alternative en- 
semble average of the scaling law, and how it performs 
concerning the symmetry relation. This latter method 
manifests itself as the most adequate technique to carry 



out the numerical MFA, leading to an even better agree- 
ment with ll]). 



II. MFA USING ENSEMBLE AVERAGE 

The numerical MFA is based on an averaged form of 
the scaling law for the giPR in the limit A = l/L -^ 0, 
where the contributions from all finite-size critical wave- 
functions are propertly taken into account. The scaling 
law for the ensemble average involves the arithmetic av- 
erage of Pq over all realizations of disorder, 



(P,(A)) ex A^° 



= (<?) 



(3) 



where (• • • ) denotes the arithmetic average over all states. 
Thus the definition of the scaling exponents is 



'{q) 



lim 



InA 



(4) 



and the corresponding definitions of a and /(a) can be 
written in a compact form as 



TABLE I: Linear system sizes L, volume V , number of eigen- 
states and total wavefunction values tpi, used for the numeri- 
cal MFA. 



L V — L samples 



t. 



20 



X 10^ 24 995 2 x 10** 



30 9 X 10'' 25 025 6.8 x 10** 

40 6.4 X lO"* 25 025 1.6 x lO'' 

50 1.3 X 10^ 25 030 3.1 x 10^ 

60 2.2 X 10^ 25 030 5.4 x 10^ 

70 3.4 X 10^ 24 950 8.6 x 10^ 

80 5.1 X 10^ 25 003 1.3 x 10^" 

90 7.3 X 10^ 25 005 1.8 x 10^" 



100 1 X 10® 


25 030 


2.5 X 10^" 


140 2.7 X 10*^ 


105 


2.9 X 10** 


160 4.1 X lO'' 


125 


5.1 X 10* 


180 5.8 X 10® 


100 


5.8 X 10** 


200 8 X 10® 


100 


8 X 10** 


210 9.3 X 10® 


105 


9.7 X 10** 


240 1.4 X 10^ 


95 


1.3 X 10^ 



|im A./^4(g,A)ln4(l,A) 



(5a) 



1 / ^" ~ 
/°- ^ /(a°-) = lim — ( Y, Sk{q, A) ln4(<Z, A) 



\fe=i 



(5b) 



Here (5/c(g, A) = /^'(A)/(Pg(A)), which is not normal- 
ized for every wavefunction but after the average over all 
of them. Let us emphasize that although Eqs. ([5]) are 
handy analytically, it is much more useful for numeri- 
cal purposes to develop them in longer expressions with 
simpler factors (see Sec.lV|. 

In contradistinction to the typical average^ which is 
determined by the behaviour of representative wavcfunc- 
tions, the ensemble average weighs the contribution of 
all wavefunctions equally, including rare (and hence not 
representative) realizations of the disorder. These rare 
events are indeed responsible for the negative values of 
/(a). Therefore it is very important to take them into 
account by doing the ensemble average, if one wants to 
have a complete picture of the singularity spectrum. We 
emphasize that in the thermodynamic limit both averag- 
ing processes must provide the same singularity spectrum 
in the positive region. The relation between typical and 
ensemble averaging has been previously commented in 
the Hteraturei^i^ 



III. DATA FOR THE 3D ANDERSON MODEL 
AT CRITICALITY 

The standard tight-binding Anderson model^^ with 
uniform diagonal disorder of mean and width Wc and 



nearest-neighbour hopping is considered in a cubic lat- 
tice of volume L^. The band width is 6 at zero disorder. 
We use the critical disorder at Wc = 16.5 and for every 
disorder realization the L^ x L^ Hamiltonian is diagonal- 
ized with periodic boundary conditions to obtain the five 
eigenstates closest to the band centre i? = 0.— The whole 
set of data used for the analysis, including system sizes 
and number of samples for each, is described in detail in 
Table [H We refer the reader to Ref. [20 for more technical 
details and to Refs. [23 and l25| for recent comprehensive 
reviews of the subject. 



IV. SCALING WITH BOX SIZE 



The easiest way to approach the thermodynamic limit 
in the scaling law ([3]) is considering the limit Z — > for the 
box size I. Using this method, we only need realizations 
for a system with a fixed linear size L, that is partitioned 
equally into an integer number of smaller boxes of linear 
size I. This is the same partitioning scheme that we have 
previously considered when studying the typical average 
of the scaling law.— 

For each state, the g-th moment of the box probability 
/i'(^) is evaluated in each box, and Pq is obtained by 
summing the contribution from all boxes. The scaling 
behaviour ([3|) is then obtained for different values of I. 
In all the computations the values of the box size ranges 
in the interval 10 ^ / ^ L/2. 
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FIG. 2: Singularity spectrum (black line) obtained using box- 
size scaling of the ensemble average of Pq for system size 
L = 100 with 2.5 x 10* states. The error bars are equal 
to one standard deviation. The corresponding symmetry- 
transformed spectrum /(2d — a) = f{a) + d — a is shown 
in black dashed line. The values for the linear correlation co- 
efficient r^ and quality-of-fit parameter Q for both a™^ and 
/°"''(a) are shown in the bottom shaded panel. 



A. General features of /''"'* (a) 




FIG. 3: (a) Mass exponents r™''(g) and (b) generalized fractal 
dimensions Dq"° corresponding to the singularity spectrum 
in Fig. [2] Dashed lines in upper panels highlight the values 
Do = d and to = —d. Symbols highlight integer values of q. 
Panel (c) : linear fits of Eq. Q . Only fits for integer values of q 
ranging from q = —5 (top) to g = 9 (bottom) are shown. The 
value of T'"^^{q) is given by the slope of the fits. Data points 
for g 7^ have been properly shifted vertically for optimal 
visualization. Data for g = highlighted with filled symbols. 
Standard deviations are contained within symbol size in all 
panels. 



The singularity spectrum for L = 100 having 2.5 x 10^ 
states is shown in Fig. [2] with its symmetry transformed 
spectrum. The first thing to notice is that the /™'^(a) 
spectrum attains negative values in the region of small a, 
corresponding to high values of the wave function ampli- 
tudes. The negative region of the multifractal spectrum 
describes the scaling of certain sets of unusual values of 
{ipil which only occur for rare critical functions. Let us 
recall that /(a) < is the fractal dimension of the set of 
points where ji/)!! r^ L~°', which implies that the number 
of such points decreases with the system size as L^'-^'"^'. 
These negative dimensions are then determined by events 
whose probability of occurrence decreases with the sys- 
tem size. The negative part of the spectrum provides 
valuable information about the distribution of wavefunc- 
tion values for a finite-size system near the critical point 
and is needed to give a complete characterization of the 
multifractal nature of the critical states at the metal- 
insulator transition. At the left-half part of /™''(a) in 
Fig. [21 we observe its termination in the negative region 
towards a — > 0. The values of a and f{a) are obtained 
from the slopes of the linear fit of Eqs. (O via a gen- 
eral x^ minimization taking into account the statistical 
uncertainty of the averaged right-hand side terms. The 
behaviour of the linear correlation coefficient r^ and the 
quality-of-fit parameter Q for the different parts of the 
spectrum (corresponding to different values of the mo- 
ments q) is shown in the bottom panel of Fig. [51 The r^ 
value is very near to one for almost all a which shows the 
near perfect linear behaviour of the data points. The pa- 



rameter Q gives an estimation on how reliable the fits are 
according to the error bars of the points involved in the 
fits. The unusual decrease of Q observed around a = 3, 
corresponding to g ~ 0.5, in Fig. [21 is due to an underesti- 
mation of the standard deviations of the points in the fits, 
since the linear correlation coefficient is still very high in 
this region. It can also be seen that the uncertainties 
for the points of /™'' (a) tend to grow when approaching 
the ends of the spectrum. This effect is more significant 
when doing ensemble average, but it should be naturally 
expected since the higher the value of \q\ is, the more the 
numerical inaccuracies of iV'iP a-re enhanced, specially 
in the region of negative q, corresponding to the right 
branch of the spectrum. The mass exponents r°"'^(g) 
and the fits of Eq. ([4]) are shown in Fig. [3l along with 
the generalized fractal dimensions D^'^^ = T"^^{q)/{q— 1) 
corresponding to the spectrum in Fig. [21 



B. Effects of system size and disorder realizations 
on /'''"'{a) 

In Fig. m we study the effects of the number of states 
and disorder realizations on /™''(q;) for L = 60 having 
2.5 X 10^ and 2.5 x 10^ states. Considering two particu- 
lar q values at each tail, when the number of samples is 
increased we see that the domain of /°"^(q;) is enlarged. 
The point corresponding to a given q appears later in the 
spectrum and thus the left end reaches more negative val- 




N 


\ 


' (b) 
<7 = -1.0 


q=-l.5 

1 


\- 



-2 



5.5 



6.5 



FIG. 4: Left (a) and right (b) branches of the singularity 
spectrum obtained using box-size scaling and ensemble aver- 
age for system size L — 60 with 2.5 x 10^ (grey) and 2.5 x 10^ 
(black) number of states. The filled symbols denote q = 5.0 
(a) and q — —1.0 (b). The empty symbols mark q = 7.0 (a) 
and q — —1.5 (b). The error bars are equal to one standard 
deviation. 
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FIG. 5: Left (a) and right (b) branches of the singularity 
spectrum obtained using box-size scaling and ensemble aver- 
age for system sizes L = 60 (grey) and L = 100 (black) where 
each has 2.5 x 10* states. The filled symbols denote q = 6.0 
(a) and q = —1.0 (b). The empty symbols mark q = 10.0 
(a) and q = —1.5 (b). The error bars indicate one standard 
deviation. 
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FIG. 6: Measure of degree of symmetry of the multifractal 
spectrum obtained from ensemble average doing scaling with 
box size. The upper panel shows the numerical evaluation of 
the symmetry law as a function of q for system sizes L = 240 
(dashed black) with 95 states, L = 60 (gray) with 2.5 x 10* 
states and L = 100 (sohd black) for 2.5 x 10* states. The bot- 
tom panel shows 5f{a) versus a. There's no correspondence 
between the abscissa axes of the upper and lower plots. For 
clarity, two values of q for the black line are explicitly written. 



number of states. Once again, we take two particular q 
values at each tail as shown in panels (a) and (b) and 
observe how their locations change when the system size 
is varied. When we consider a bigger system size with 
the same number of reahzations, the domain of /''"''(a) 
tends to decrease, and so for the same q range, we see 
less negative values at the left end [Fig. [5ja)]. In other 
words to be able to observe the same extent of the neg- 
ative f°'^^{a) values of L = 60, one must average over 
more states when a bigger system size such as L = 100 
is considered. The same shrinking effect also occurs in 
the right branch of the spectrum. This unexpected be- 
haviour is due to the nature of the ensemble averaging 
process, that is strongly determined by the contribution 
of rare events which are less likely to happen for larger 
systems. This important effect will be discussed in more 
detail in Sec. IVl 



C. Symmetry relation 



ues with more states [Fig.Ul^a)]. The same stretching ef- 
fect can also be observed for the right branch in Fig.[4ljb). 
Additionally the reliability of the data points in the sin- 
gularity spectrum is significantly improved as shown by 
the huge decrease in their uncertainties. These effects 
prove the strong dependence of the ensemble averaging 
on the number of samples taken. 

The effect of the system size on the shape of the singu- 
larity spectrum is presented in Fig. \E\ Here we consider 
system sizes L = 60 and L = 100 each having 2.5 x 10* 



In the upper panel of Fig. [6] we give a numerical evalua- 
tion of the symmetry law. An approximate estimation of 
the symmetry law is also shown in the lower panel using 
([2]) ^ which measures the distance between the spectrum 
and its symmetry-transformed counterpart. We compare 
data fori = 240 (95 states), L = 60 (2.5x10* states) and 
L = 100 (2.5 X 10* states). Our results show that in gen- 
eral the closest agreement to the symmetry in the singu- 
larity spectrum is achieved for the cases with the highest 
number of disorder realizations, in particular for L = 100 
(/(a) shown in Fig. [5]). Although around the symmetry 



point q — 1/2 the spectrum obtained using the largest 
system size available, L = 240 with 95 states, tends to 
behave slightly better (inset in upper panel of Fig. [2]), the 
tendency is inverted when looking at a broader region of 
q values. This result is a clear manifestation of how im- 
portant the number of disorder realizations is when doing 
ensemble average. 

It is therefore clear that the more realizations, the bet- 
ter the symmetry is. Obviously a bigger system size helps 
reduce finite-size effects, but we have shown that increas- 
ing system size also implies generating more states in or- 
der to obtain the same extent of /™''(a). Thus from a 
numerical viewpoint an agreement between system size 
and disorder realizations must be found to optimize the 
use of box-size scaling and ensemble averaging. 
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V. SCALING WITH SYSTEM SIZE 

The scaling with the system size may be the most ad- 
equate way to describe the thermodynamic limit of the 
scaling law for the gIPR ([3]) (L — > cx)), however the nu- 
merical eigcnstate problem is highly demanding for very 
large 3D systems i^ 

The formulae @ and ([5]) for the singularity spec- 
trum are now affected by the substitution: limA_»o =^ 
— liniL^oo- As for the typical average^^ the box size / 
which determines the integrated probability distribution 
is set to ^ = 1 for non-negative moments {q ^ 0) and to 
a value I > 1 (usually I = 5) for g < 0, in order to mini- 
mize the errors and the uncertainties in the right branch 
of /(a). For the case I = 1 the formulae to obtain the 
spectrum reduce to 



-a™"lnL^ 



-f 



^'^InL 



(EJV'.P^ln 


IV'.P) 







(6a) 



EJ^.P^ 



-In (El'/'' 



2q 



(6b) 



We note the clear difference between the ensemble av- 
erage ( El) a nd the typical average techniques [Eqs. (14) 
in RefH^I. The values of a™'' and /™" are obtained 
from the slopes of the linear fits of the averaged terms in 
Eqs. ® (q ^ 0) and © (g < 0) versus luL, for different 
values of the system size L. 



A. General features of ^"^(a) 

The multifractal spectrum obtained from the ensemble 
average is shown in Fig. [71 where we have considered 9 
different linear system sizes ranging from L = 20 to 100 
for the scaling after averaging over ~ 2.5 x lO"' states 
for each size as shown in Table [H The branch of negative 
values characterizing /°"*'(a) can be clearly seen. The ab- 
sence of an infinite slope in the spectrum when crossing 



FIG. 7: Singularity spectrum obtained from ensemble aver- 
aging. System sizes from i = 20 to 100 have been consid- 
ered with ~ 2.5 X 10* different wavefunctions for each system 
size as shown in Table [H The dashed line is the symmetry- 
transformed spectrum. The values of q range from q = —2 to 
g = 7 with a step of 0.1 (^ = 1 for g ^ and / = 5 for g < 0). 
Symbols highlight the values corresponding to integer q. Er- 
ror bars in grey are standard deviations. The lower panel 
shows the linear correlation coefficient (r^) and the quality- 
of-fit parameter (Q) of the hnear fits to obtain the values for 
a and /(a). 



the abscissa axis must also be emphasized. As discussed 
in Ref. l20J this confirms the fact that the divergence of 
the slope at the termination points observed when do- 
ing the typical average, /*''P(a), is purely a finite-size ef- 
fect, since both averages must provide the same result for 
/(a) ^ in the thermodynamic limit. This is also sup- 
ported by the systematic shift of the left end of /*^p(q;) 
to smaller values of a whenever more states or larger sys- 
tem sizes are considered)^ The error bars for the values 
of /°"''(a) are considerably larger than the ones obtained 
for /'''p(q;) using the same system sizes and disorder re- 
alizations [Fig. (7) in Ref. |20|. This is of course a conse- 
quence of having larger errors for the points used in the 
linear fits, shown in Fig. [8] These higher uncertainties are 
due to the nature of the average itself and the probability 
distribution function for the generalized IPR. The proba- 
bility density for Pq is an asymmetric function around its 
maximum with long tails j^ii^ resembling a log-normal 
distribution. The calculation of the arithmetic average 
of Pq , involved in the ensemble average is therefore much 
more heavily based on the number of disorder realiza- 
tions than the determination of the geometric mean used 
for the typical average, and thus larger uncertainties and 
slower convergence must be expected for the ensemble- 
averaged situation with the same number of wavefunc- 
tions. Regarding the errors in the values of /°"*'(a), it 
is remarkable how their magnitude grows, for high val- 
ues of |g|, apparently at the same rate as the spectrum 
deviates from the symmetry-transformed curve (dashed 
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FIG. 8: Linear fits of Eqs. dS} for a°"" values (left) and f, 
values (right) of the singularity spectrum in Fig. [T] Only fits 
for integers values of q ranging from g = 7 (top) to g = — 2 
(bottom) are shown. The values of q^"^ and /q"° are given by 
the slopes of the fits. Data points for q j^ have been properly 
shifted vertically to ensure optimal visualization. Data for q — 
highlighted with filled symbols. When not shown, standard 
deviations are contained within symbol size. 



line in Fig. [7]). This suggests that it might be possible to 
observe almost a perfect agreement with the symmetry 
law, using this small range of system sizes for the scaling, 
if the number of realizations were large enough. 



B. Effects of the number of disorder realizations on 

reus I \ 

f («) 

The effect of increasing the number of states in the 
ensemble average can be seen in Fig. ^ for scaling with 
L S [20,100]. A reduction of the standard deviations 
must be expected whenever more realizations are taken 
into account. To make this clear we have considered two 
situations: averaging over 10'^ states or over ~ 2.5 x lO'' 
states for each size. In Fig. [9] the points with the speci- 
fied vertical uncertainty appear later (for higher values of 
q) on the left and right branches of the spectrum when 
we increase the number of states in the average. This 
clearly means that for a fixed position on the f°'^^{a) 
curve, the uncertainty shrinks when more states are in- 
cluded. There is however, another significant effect that 
must be emphasized. In Fig. [9l^a)(b), the spectrum ob- 
tained for values of q higher than the ones indicated is 
represented by dashed lines. For the average including 
only 10'^ samples for each size, the values of the spec- 
trum for high l^l are completely absurd and f{a) behaves 
in an unexpected way, showing "kinks" and bumps as a 
consequence of a loss of precision in the fits caused by 
very large uncertainties. This implies that by increasing 
the number of states in the ensemble average not only 
the standard deviations are reduced for each point, but 
the domain of accessible values for /(a) is also enlarged, 
e.g. for more wavefunctions the spectrum reaches more 



FIG. 9: Left (a) and right(b) branches of the singularity spec- 
trum obtained from ensemble averaging scaling with system 
sizes from L = 20 to 100 and same number of states for each 
size: (grey) 10"^, (black) ~ 2.5 x 10*. The values of q range 
from g = — 10 to g = 10 with a step of 0.1 (/ = 1 for g ^ and 
I — 5 for I < 0). The vertical standard deviation for the points 
marked with filled symbols is always cry = 1.0, and only the 
uncertainty for a has been included, for clarity. The q value 
corresponding to each of the symbols is indicated. Dashed 
lines represent the spectrum in each case for higher values of 
q. Inset (c) shows the change in the value /(a = 6) for the 
cases: (grey) 5 x 10^ states for each size, (black) ~ 2.5 x 10* 
states for each size. 



negative values [Fig. [DJa)]. This is in marked contrast 
to the typical average»2£ for which one obtains almost 
the same range of f{a) independently of the number of 
states considered. When doing the ensemble average, the 
appearance of these "kinks" in the spectrum, cither for 
system-size or box-size scaling where they have also been 
observed, is always the fingerprint of a lack of sampling 
of the distributions, i.e. not enough disorder realizations. 
In the inset (c) of Fig. [21 we have also illustrated the 
behaviour of the value /(a = 6), at the upper boundary 
required by the symmetry relation, when the number of 
states is changed from 5 x lO'^ to ^ 2.5 x 10* for each 
system size. The spectrum tends to be in better agree- 
ment with the upper boundary required by ^ when the 
number of disorder realizations increases. 



C. Effects of the range of system sizes on f'^'^^{a) 

In order to study the effects of the system size, we show 
in Fig. [TU] the singularity spectrum obtained doing scal- 
ing in different intervals: L e [40, 100] and L e [140, 240], 
taking ~ 100 wave functions for each size in both cases. 
The fact that we have only averaged over 100 states for 
each size, makes the standard deviations noticeably large, 
however this does not affect the conclusions qualitatively. 
When we consider larger system sizes for a fixed number 
of disorder realizations, the region where we can reliably 
obtain the multifractal spectrum shrinks. Moreover if we 
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FIG. 10: Singularity spectrum obtained from ensemble aver- 
aging scaling with 7 system sizes L £ [40, 100] and 10^ states 
for each size (grey), and 6 system sizes L G [140, 240] and 
~ 10 states for each size (black). The values of q range from 
q = — 10 to g = 10 with a step of 0.1 (Z = 1 for q ^ Q and 
I = 5 for I <Q). Filled symbols correspond to points with the 
same vertical standard deviation (cr/ ~ 1.0, not shown for 
clarity) . Dashed lines represent the spectrum in each case for 
values of q higher than the ones indicated beside their corre- 
sponding points. Inset (a) shows the spectrum for a different 
set of data: L £ [20, 60] (grey) and L £ [60, 100] (black) with 
~ 2.5 X lO'* states for each size in both cases. 
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FIG. 11: Measure of degree of symmetry of the multifractal 
spectrum obtained from ensemble average doing scaling with 
system size. The upper panel shows the numerical evaluation 
of the symmetry law as a function of q. Dashed black: 7 
system sizes from L = 40 to 100 and 10^ states for each. 
Grey: 6 system size from L — 140 to 240 and ~ 10^ states 
for each. Solid black: 9 system sizes from L = 20 to 100 and 
~ 2.5 X 10* states for each. The bottom panel shows Sf(a) 
versus a. There's no correspondence between the abscissa 
axes of the upper and lower plots. For clarity, two values of 
q for the black line are explicitly written. 



go to high enough values of jql (highlighted by dashed 
lines in Fig. [TU)) . it can be noticed how the wrong be- 
haviour of /(a) is enhanced. This is a very counterin- 
tuitive result, since one would expect that for increasing 
system sizes, the number of disorder realizations needed 
to obtain the spectrum with a given degree of reliability 
should decrease proportionally - that is in fact what hap- 
pens with the typical averaging.— However for ensemble 
averaging the conclusion is just the opposite: if you want 
to improve the spectrum in a given region of the tails 
and you consider larger system sizes to reduce finite size 
effects, the number of disorder realizations must also be 
increased. This is due again to the nature of the ensem- 
ble averaging process and the shape of the distributions 
for the gIPR. The arithmetic average is heavily based on 
rare events which are less likely to appear for bigger sys- 
tems, and so the number of realizations has to grow with 
the system size in order to include the proper amount 
of rare events. This can be more clearly understood in 
the region of negative fractal dimensions. We know that 
the number of points in a single wave function such that 
lipfl ~ L-" where f{a) < 0, is L-l-'^^")! < 1. Therefore 
to be able to see the region of negative fractal dimensions 
we would need a number of states Af such that we can 
guarantee TVi^i-''-"" ^ i This implies that the number 
of disorder realizations must go as A/" ~ il-^^")', and thus 
it increases with the system size. This effect can be ob- 
served in the inset (a) of Fig. 1101 where we have compared 
scaling with sizes L e [20, 60] and L S [60, 100] with 
~ 2.5 X lO'' states for each size in both cases. For higher 



sizes and the same number of states, we are not able to 
see the same region of negative fractal dimensions. Aside 
from this effect, it must nevertheless be emphasized that 
when we consider larger system sizes, the right branch of 
the spectrum tends to find a better agreement with the 
upper boundary required by the symmetry law. 



D. Symmetry relation 

To discuss the fulfilment of the symmetry using en- 
semble average and scaling with system size, let us look 
at Fig. I 111 where the numerical evaluation of the symme- 
try law ll]) is shown for different ranges of system sizes 
and disorder realizations. The best result, according to 
the symmetry, corresponds undoubtedly to the case with 
the highest number of disorder realizations, ~ 2.5 x 10^, 
for which the scaling analysis involves sizes from L ~ 20 
to 100. The difference is remarkable between the situ- 
ation corresponding to (i) L E [40, 100] averaging over 
100 states only, where the symmetry is hardly satisfied 
at all, and (ii) the best case where the development of the 
plateau for aq+ai-q^2d around q = 0.5 can be seen very 
clearly. The spectrum obtained for L £ [120, 240] with 
~ 100 states for each size also deviates noticeably from 
the symmetry. These differences can also be seen in the 
bottom panel of Fig. I 111 where the degree of symmetry is 
estimated by 6f(a). 

For the ensemble average going to very large system 
sizes is not the best strategy unless one can generate 



an increasing number of states. It must be clear that 
of course finite size effects will be reduced using large 
sizes but the number of wavefunctions used for the av- 
erage has to grow with the system size considered. For 
a given range of system sizes, increasing the number of 
states improves the reliability of data, enlarges the acces- 
sible domain of /(a), specially in the region of negative 
dimensions, and improves the symmetry. 



VI. COMPARISON OF DIFFERENT SCALING 
AND AVERAGING APPROACHES 

Taking the symmetry relation ([1]) as a measure of the 
quality of the numerical MFA, let us compare the re- 
sults of the different scaling and averaging techniques. In 
Fig- EH we show the best analyses obtained from box-size 
scaling and system-size scaling using typical average and 
ensemble average in both cases. Data corresponding to 
the typical average has been extracted from Ref. [20. The 
performance of the system-size scaling technique (solid 
lines) is clearly much better than box-size scaling (dashed 
lines). This is not a very surprising result, since one ex- 
pects finite-size effects to be more enhanced in box-size 
scaling. For each of the scaling procedures the ensemble 
average (black) is also better than the typical average 
(grey). This may be not be so intuitive, since due to the 
nature of the distribution functions for Pg ^L^ one might 
expect the typical average to be a better choice. How- 
ever it turns out that the ensemble average does better in 
revealing the true behaviour in the thermodynamic limit. 

Let us recall the strategies that give the best result 
for each of the techniques. For typical average the best 
symmetry is achieved using the largest system sizes avail- 
able, either for box-size or system-size scaling, although 
the number of realizations is not the highest 3^ On the 
other hand, using ensemble average, the safest choice is 
to consider smaller system sizes for which a very large 
number of disorder realizations can be obtained. 



VII. CONCLUSIONS AND OUTLOOK 

In this work we have studied the symmetry law ([1]) 
for the multifractal spectrum of the electronic states at 
the metal-insulator transition in the 3D Anderson model, 
using the ensemble-averaged scaling law of the gIPR ([3]). 
A detailed analysis has revealed how the MFA is affected 
by system size and number of samples. System-size scal- 
ing with ensemble average has manifested itself as the 
most adequate method to perform numerical MFA, in 
contrast to box-size scaling and typical average which had 
been mainly the method of choice in previous studies.— i^i^ 
Since the ensemble average is strongly based on the num- 
ber of disorder realizations, from a numerical point of 
view, the best strategy to carry out the analysis is to 
consider a sensible range of system sizes for which a very 
large number of states can be generated. 
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FIG. 12: Comparison of degree of symmetry for spectra 
obtained using: box-size scaling typical average [BS-TYP] 
(dashed grey), box-size scaling ensemble average [BS-ENS] 
(dashed black), system-size scaling typical average [SS-TYP] 
(solid grey), system-size scaling ensemble average [SS-ENS] 
(solid black). The best spectrum for each case has been con- 
sidered: [BS-TYP] L = 240 (95 states), [BS-ENS], L = 100 
(2.5 X 10* states), [SS-TYP] L € [140,240] (- 10^ states for 
each size), [SS-ENS] L G [20, 100] (~ 2.5 x lO" states for each 
size) . Data for typical average extracted from Ref. I20I . The 
upper panel shows the numerical evaluation of the symmetry 
law as a function of q. The bottom panel shows Sf{a) versus 
a. There's no correspondence between the abscissa axes of 
the upper and lower plots. 



All our results suggest that the symmetry law is true 
in the thermodynamic limit, since a better agreement is 
found whenever a high enough number of disorder re- 
alizations and larger system sizes are considered. The 
symmetry relation ^ then provides a powerful tool to 
obtain a complete picture of /(a) at criticality, since it 
would suffice to obtain numerically the spectrum in the 
most reliable region, g > 0, and apply the symmetry to 
complete the function for q < 0. 

The results obtained for f{a) also provide some use- 
ful information about the validity of previous analytical 
results. The perturbative analysis in d = 2 + e made by 
Wegner— led to the following spectrum 



,r{a)^d- 



[a^id + ef 



C(3) 
64 



4e 
[{a - df 



(7) 



[a — dy 



15e^ 



where (^{x) denotes the Ricmann zeta function. The first 
two terms in ^ constitute the usual parabolic approxi- 
mation (PA). The extra quartic term is an overestimation 
in 3D (e = 1), as explicitly stated by Wegner, which gives 
a non-acceptable spectrum as can be seen in Fig. [T31 To 
obtain the correct /(a) at e = 1 all the other terms in 
the perturbation series are required. Remarkably Wcg- 
ner's result obeys the symmetry relation ([5|), provided 
the spectrum is indeed terminated at a = and a = 2d. 
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FIG. 13: Singularity spectra given by the Parabolic Approxi- 
mation (dashed gray), Eq. (dashed black), numerical cal- 
culation doing scaling with system-size and ensemble average 
described in Fig. [7] (solid gray) and conjectured /(q) at the 
MIT for the 3D Anderson model (solid black). The dotted 
black areas highlight forbidden regions for /(a). The grey 
shaded and white areas can be connected through the sym- 
metry relation (HI,©- 



The symmetry can also be easily checked by calculating 
Aq, which can be written as a power series in e whose 
coefficients are all invariant under the transformation 
q ^r (\ — q). The results from the numerical MFA seem 
to suggest that /(a) tends to approach a finite negative 
value when a — > 0. This would imply that the spectrum 
has a termination point. The existence of a termination 
point at a = requires that the function t{(i) has a 



zero slope for q greater than a certain critical value qc^ 
This latter fact is also suggested in Fig. [3{a). Moreover 
the numerically obtained spectrum is between the PA 
and /^(a) as shown in Fig. 1131 In view of this we are 
tempted to conjecture that perhaps the PA and /^(a) 
arc bounds for the spectrum in the thermodynamic limit. 
In this were true then /(a) must necessarily terminate 
at finite values. Furthermore let us speculate that the 
left termination point is j(a = 0) = —d which implies 
/(a = 2d) = according to the symmetry [Fig. [T^ . 
and thus the spectrum would not attain negative values 
at its right end, determined by the smallest wavefunc- 
tion amplitudes. The resulting /(a) at the 3D Anderson 
transition would hence be close to the black line shown 
in Fig. [HI 

In conclusion, in spite of the large amount of informa- 
tion that the numerical analyses, here and in Ref. |20| . 
and the symmetry relation ^ have provided, the com- 
plete picture of the multifractal spectrum for the MIT 
in 3D is still elusive. In particular further research is 
needed to confirm the existence of termination points and 
whether these happen at negative values on both sides, 
since this has important implications upon the distribu- 
tion functions of the wavefunction amplitudes near the 
localisation-dclocalisation transition. 
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